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An attempt has been made to extract the evolution of radial flow from the analysis of the ex- 
perimental data on electromagnetic probes measured at SPS and RHIC energies. The transverse 
momentum (pt) spectra of photons and dileptons measured by WA98 and NA60 collaborations 
respectively at the CERN-SPS and the photon and dilepton spectra obtained by PHENIX collabo- 
ration at BNL-RHIC have been used to constrain the theoretical models. We use the ratio of photon 
to dilepton spectra to extract the flow, where some model dependence are canceled out. Within the 
ambit of the present analysis we argue that the variation of the radial velocity with invariant mass 
is indicative of a phase transition from initially produced partons to hadrons at SPS and RHIC 
energies. 
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I. INTRODUCTION 

The hot and dense matter expected to be formed 
in the partonic phase after ultra-relativistic heavy 
ion collisions dynamically evolve in space and time 
due to high internal pressure. Consequently the sys- 
tem cools and reverts back to hadronic matter from 
the partonic phase. Just after the formation, the en- 
tire energy of the system is thermal in nature and 
with progress of time some part of the thermal en- 
ergy gets converted to the collective (flow) energy. 
In other words, during the expansion stage the total 
energy of the system is shared by the thermal as well 
as the collective degrees of freedom. The evolution 
of the collectivity within the system is sensitive to 
the Equation of State (EoS). Therefore, the study 
of the collectivity in the system formed in the quark 
gluon plasma (QGP) phase will be useful to shed 
light on the EoS ^iHj] (see 0, Q for review) and on 
the nature of the transition that may take place dur- 
ing the evolution process. It is well known that the 
average magnitude of radial flow at the freeze-out 
surface can be extracted from the transverse momen- 
tum (pt) spectra of the hadrons. However, hadrons 
being strongly interacting objects can bring the in- 
formation of the state of the system when it is too 
dilute to support collectivity i.e. the parameters 
of collectivity extracted from the hadronic spectra 
are limited to the evolution stage where the collec- 
tivity ceases to exist. These collective parameters 
have hardly any information about the interior of the 
matter. On the other hand electromagnetic (EM) 
probes, ie.photons and dileptons are produced and 
emitted [l-Q (see [9l-[Tl| for review) from each space 
time points. Therefore, estimating radial flow from 
the EM probes will shed light on the time evolution 
of the collectivity in the system fl^l. 

The invariant momentum distribution of photons 
produced from a thermal source depends on the tem- 
perature (T) of the source through the thermal phase 
space distributions of the participants of the reac- 



tion that produces the photon [l3|. As a result the 
Pt spectra of photon reflects the temperature of the 
source. Hence ideally the photons with intermedi- 
ate Pt values (~ 2 — 3 GeV, depending on the value 
of initial temperature) reflect the properties of QGP 
(realized when T > Tc, Tc is the transition temper- 
ature). Therefore, one should look into the pt spec- 
tra for these values of pt for the detection of QGP. 
However, for an expanding system the situation is 
far more complex. The thermal phase space factor 
changes - several factors e.g. the transverse kick re- 
ceived by low Pt photons due to flow originating 
from the low temperature hadronic phase (realized 
when T <Tc) populates the high pT part of the spec- 
tra [IJ]. As a consequence the intermediate or the 
high Pt part of the spectra contains contributions 
from both QGP and hadrons. For dilepton the sit- 
uation is, however, different because in this case we 
have two kinematic variables - out of these two, the 
Pt spectra is affected by the flow, however, the pT in- 
tegrated invariant mass (M) spectra is unaltered by 
the flow in the system. It should be mentioned here 
that for M below p peak and above (f) peak dilep- 
tons from QGP dominates over its hadronic coun- 
terpart (assuming the contributions from hadronic 
cocktails are subtracted out) within the framework 
of the present model. However, the spectral func- 
tion of low mass vector mesons (mainly p) may shift 
toward lower invariant mass region due to non-zero 
temperature and density effects. As a consequence 
the contributions from the decays of such vector 
mesons to lepton pairs could populate the low M 
window and may dominate over the contributions 
from the QGP phase [l^ (and also El for re- 
view). In the present work such thermal effects are 
not considered. All these suggests that a judicious 
choice of pT and M windows will be very useful to 
characterize the flow in QGP and hadronic phase. 
However, there are still some difficulties. The calcu- 
lations of EM probes from thermal sources depend 
on the parameters like initial temperature (T^), ther- 
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malization time (t^), chemical freeze-out tempera- 
ture (Teh), kinetic freeze-out temperature (Tf) etc, 
which are not known uniquely. To minimize the de- 
pendence of thermal sources on these parameters the 
importance of the ratio of the transverse momentum 
spectra of photon to dilepton has been emphasized 
previously (T6l - [l8| in order to overcome the above 
mentioned uncertainties. It may be mentioned here 
that in the limit of M — > the lepton pairs (vir- 
tual photons) emerge as real photons. Therefore, the 
evaluation of the ratio of the pt spectra of photons to 
dileptons for various invariant mass bins along with 
a judicious choice of the pt and M windows will 
be very useful to extract the properties of QGP as 
well as that of hadronic phase. This will be demon- 
strated in the present work by analyzing WA98 and 
PHENIX photons and NA60 and PHENIX dilepton 
spectra. 

The paper is organized as follows. In the next sec- 
tion, the production of thermal photons and dilep- 
tons are briefly outlined. In section III the expansion 
dynamics of the system is discussed. Section IV is 
devoted to results and discussions and finally in sec- 
tion V, we present summary and conclusions. 



II. PHOTONS AND DILEPTONS 
PRODUCTION 

The Pt and M spectra of EM probes measured 
in the experiments are mingled with all the sources 
of productions - broadly categorized as: (i) prompt 
productions resulting from the interactions of the 
partons of the colliding nuclei, (ii) thermal produc- 
tions from the interactions of thermal partons as 
well as from thermal hadrons and (iii) finally from 
the decays of the long (compared to strong inter- 
action time scale) lived mesons. The contributions 
from pp collisions at a given collision energy can 
be used as a bench-mark to estimate the hard con- 
tributions. To estimate the thermal contributions 
we adopt the following procedure: thermal contribu- 
tions = contributions from heavy ion collision minus 
NcoU X contributions from pp collisions where Ncoii 
is the number of nucleon-nucleon interactions in the 
nuclear collisions at a given centrality. Instead of 
evaluating the hard contributions by applying pQCD 
we use the experimental data both for photons and 
dileptons from the pp collisions wherever available 
to minimize the uncertainties in the contributions of 
category (i). 

For collisions with large nuclei e.g. Au + Au or 
Ph -\- Pb, the number of valence d quarks are more 
than the number of valence u quark because in these 
nuclei there are more neutrons than protons. The 
magnitude of electric charge of the d quark is half of 
that of u quarks. Consequently, in the production of 



EM probes from Au + Au collisions - a larger num- 
ber d quarks with lesser charge and comparatively a 
smaller number of u quarks with larger charge are in- 
volved. Therefore, for the photon production from 
Au + Au interaction for the category (i) is not a 
mere superposition of the yield from p + p interac- 
tion. However, in the present work we concentrate 
in the kinematic region of central rapidity where the 
number of valence quarks are small. In fact, this 
is negligible for RHIC energy. It is worth mention- 
ing here that the role of EM probes produced from 
other mechanisms - like fragmentation of high energy 
quarks and the interaction of high energy partons 
with thermal QGP medium |il9i] are ignored here. 
It is expected that at the px domain of our inter- 
est omission of these mechanism will not change the 
final results significantly. 



A. Production of thermal photons and 
dileptons 

For the present work photons and dileptons from 
thermalized matter of partons and hadrons play the 
most crucial role. The rate of thermal dilepton pro- 
duction per unit space-time volume per unit four 
momentum volume is given by 
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where a is the EM coupling constant, Imll^ is the 
imaginary part of the retarded photon self energy 
and fBE(,E,T) is the thermal phase space distribu- 
tion for Bosons. L(p^)=(l -I- 2vr? IjP'')^\ — Am? jp^, 
arises from the final state leptonic current involving 
Dirac spinors of mass m. As mentioned before, in 
the limit of vanishing M a lepton pair i. e. a virtual 
photon appears as a real photon. Therefore, the 
real photon production rate can be obtained from 
the dilepton emission rate by replacing the product 
of EM vertex 7* — > the term involving final 

state leptonic current and the square of the (vir- 
tual) photon propagator by the polarization sum for 
the real photon. For an expanding system E should 
be replaced by u^p'^, where p^ and u'^ are the four 
momentum and the hydrodynamic four velocity re- 
spectively. 



B. Thermal photons 

The Hard Thermal Loop [20| approximations has 
been used by several authors \2V\ to evaluate the 
photon spectra originating from a thermal source 
of quarks and gluons. The complete calculation of 
emission rate of photons from QGP to order 0(aQ;s) 
has been done by resuming ladder diagrams in the 
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effective theory [22|, which has been used in the 
present work. A set of hadronic reactions with ah 
possible isospin combinations have been considered 
for the production of photons [23l - [25j from hadronic 
matter. The effect of hadronic dipole form factors 
has been taken into account in the present work as 
in (H. 



C. Thermal dileptons 

The lowest order process producing lepton pairs is 
q and q annihilation. The correction of order agO^ 
to the lowest order rate of dilepton production from 
QGP has been calculated in [2y, [23|, which is con- 
sidered in the present work. For the low M dilepton 
production from the hadronic phase the decays of 
the light vector mesons p, u and (j) has been consid- 
ered in [m HI] • The continuum part of the vec- 
tor mesons spectral functions constrained by exper- 
imental data have been included here. As men- 
tioned before the contributions from the QGP phase 
dominates the M spectra of the lepton pairs below 
p-peak and above the 0-peak if no thermal effects 
of the spectral functions of the vector mesons (see 
[l3,liil,[2i] for review) are considered. Since the con- 
tinuum part of the vector meson spectral functions 
are included in the current work the processes like 
four pions annihilations [30j are excluded to avoid 
double counting. 
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FIG. 1: Transverse momentum spectra of photon at SPS 
energy for Pb+Pb collision at mid rapidity. 
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FIG. 2: Transverse momentum spectra of photons at 
RHIC energy for Au-Au coUision for different centralities 
at mid-rapidity. 



III. EXPANSION DYNAMICS 

The space time evolution of the system formed 
in heavy ion collisions has been studied by using 
relativistic hydrodynamics with longitudinal boost 
invariance [3lj and cylindrical symmetry fs^. We 
assume that the system reaches equilibration at a 
time Ti after the collision. The initial temperature, 
Ti can be related to the measured hadronic multi- 
plicity (dN/dy) by the following relation: 



1 



1 dN 



45C(3) 4ae// nRj, dy 
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where Ra is the radius of the system, ^(3) is the 
Riemann zeta function and Oe// = 7r^(7e///90 where 
g^jf (=2x8-h7x2x2x3x Np/S) is the degen- 
eracy of quarks and gluons in QGP, 7Vi?=number of 
flavours. The value of dN / dy for various beam ener- 
gies and centralities are calculated from the following 
equation [33| : 
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and tabulated in table I. Ncoii is the number of 
collisions and contribute to x fraction to the mul- 
tiplicity dupp/dy measured in pp collision. The 
number of participants, Npart contributes fraction 



part 

The values of Npart and Ncoii 



(1 — x) of dupp/dy. 
are estimated by using Glauber Model and the re- 
sults are in agreement with [s^. We have used 
dupp/dy = 2 A3 and x = 0.1 at ^/snn = 200 GeV. It 
should be mentioned here that the values of dN/dy 
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TABLE I: The values of various parameters - thermal- 
ization time (r^), initial temperature (Ti) and hadronic 
multiplicity dN/dy - used in the present calculations. 



\/snn 


centrality 


dN 
dy 


Tiifm 


) r,(MeV) 


17.3 GeV 


0-06% 


700 


1.0 


200 


200 GeV 


0-20% 


496 


0.6 


227 




20-40% 


226 


0.6 


203 




min. bias 


184 


0.6 


200 



(through Npart and Ncoii in Eq. [3]) and hence the 
Ti (through dN/dy in Eq. [2|) depend on the cen- 
trality of the collisions. The values of Ra for dif- 
ferent centralities have been evaluated by using the 
equation Ra ~ l.liVp^^j. Some comments are in or- 
der here regarding the flow for non-central collisions 
vifhich results in non-cylindrical geometrical shape 
of the system formed after the collisions. In the 
present work the non-centrality of the nuclear col- 
lisions is reflected in the initial temperature through 
hadronic multiplicity. The geometry of the system 
due to non-centrality, should in principle be treated 
by (2-1-1) dimensional fss'l or for more rigorous treat- 
ment (3-t-l) dimensional hydrodynamical [s^ evolu- 
tion. However, we expect that the results obtained 
in the present work will not be affected substantially 
due to non-cylindrical geometric shape of the sys- 
tem. Because the flow has been extracted from the 
ratio of photon to dilepton spectra for a given cen- 
trality. 



We use the EoS obtained from the lattice QCD 
calculations by the MILC collaboration 37]. We 
consider kinetic freeze out temperature, Tf=\AQ 
MeV for all the hadrons. The ratios of various 
hadrons measured experimentally at different ^snn 
indicate that the system formed in heavy ion col- 
lisions chemically decouple at Teh which is higher 
than Tf which can be determined by the transverse 
spectra of hadrons [s^. Therefore, the system re- 
mains out of chemical equilibrium from Teh to Tf. 
The deviation of the system from the chemical equi- 
librium is taken in to account by introducing chem- 
ical potential for each hadronic species. The chem- 
ical non-equilibration affects the yields through the 
phase space factors of the hadrons which in turn af- 
fects the productions of the EM probes. The value of 
the chemical potential has been taken in to account 
following Ref. Q. It is expected that the chemical 
potentials do not change much for the inclusion of 
resonances above A. 



IV. RESULTS AND DISCUSSION 

A. pt distributions of photons and dileptons 

The prompt photons and dileptons (Drell-Yan) 
are normally estimated by using perturbative QCD. 
However, to minimize the theoretical model depen- 
dence here we use the available experimental data 
from pp collisions to estimate the hard photon and 
dilepton contributions in heavy ion collisions. The 
WA98 photon spectra from Pb-|-Pb collisions is mea- 
sured at y/SNN = 17.3 GeV. However, no data 
at this collision energy is available for pp interac- 
tions. Therefore, prompt photons for p+P collision 
at ^/SNN — 19.4 GeV has been used 



[39| to esti- 
mate the hard contributions for nuclear collisions at 
y/SNN = 17.3 GeV. Appropriate scaling (see [13] 
for details) has been used to obtain the results at 
^/sWn = 17.3 GeV. For the Pb-|-Pb collisions the 
result has been appropriately scaled by the number 
of collisions at this energy (this is shown in Fig. [T] 
as prompt photons) . The high px part of the WA98 
data is reproduced by the prompt contributions rea- 
sonably well. At low px the hard contributions un- 
der estimate the data indicating the presence of a 
thermal source. The thermal photons with initial 
temperature = 200 MeV along with the prompt 
contributions explain the WA98 data well (Fig. [ij , 
with the inclusion of non-zero chemical potentials 
for all hadronic species considered H^see also [4l|). 
In some of the previous works |42h47| the effect of 
chemical freeze-out is ignored. As a result either 
a higher value of Ti or a substantial reduction of 
hadronic masses in the medium was required p^ . 
In the present work, the data has been reproduced 
without any such effects. 

Following a similar procedure, the data [4^ from 
Au-Au collisions at ^/JWn = 200 GeV has been re- 
produced well by adding the prompt contributions 
(which is constrained by pp data at the same energy) 
to the thermal photons. The reproduction of data is 
satisfactory (Fig. [5]) for all the centralities with the 
initial temperature shown in table I (see also [49j . 

The transverse mass distribution of dimuons pro- 
duced in In-|-In collisions at ^^snn = 17.3 GeV has 
been evaluated for different invariant mass ranges 
(see [50| for details). The quantity dN/MrdMr has 
been obtained by integrating the production rates 
over invariant mass windows M^in to M^ax and 
Mt is defined as y/ < M +p'^ where < M >= 
{M^in + M^ax)/"^- The results are compared with 
the data obtained by NA60 collaborations [5l| at 
SPS energy (Fig. [3]). Theoretical results contain con- 
tributions from the thermal decays of light vector 
mesons (p, uj and 0) and also from the decays of 
vector mesons at the freeze-out 52| of the system 
has also been considered. The non-monotonic varia- 
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FIG. 3: Transverse mass spectra of dimuons in In+In 
collisions at SPS energy. Solid lines denote the theoreti- 
cal results. 



PHENIX collaboration [53| could be reproduced well 
with the same initial condition that reproduces pho- 
ton spectra [i^. In fact, the reproduction of data for 
the higher mass windows 0.5 <M(GeV)< 0.75 and 
0.81 <M(GeV)< 0.99 do not need any normalization 
factors (Fig. |4]). For lower mass windows slopes are 
reproduced well but fail to reproduce the absolute 
normalization. Therefore, it should be clarified here 
that the theoretical results shown in Fig. U for lower 
mass windows (to be precise for 0.1 <M(GeV)< 0.2, 
0.2 <M(GeV)< 0.3 and 0.3 <M(GeV)< 0.5) con- 
tain arbitrary normalization constant. However, it 
should also be mentioned at this point that for the 
extraction of the flow within the present approach 
the absolute normalization is not essential, which 
is essential is the slope (Eq. [5]). Therefore, the non- 
reproduction of the absolute normalization of the pt 
spectra of lepton pairs for the lower mass windows 
does not affect the extraction of the magnitude of 
the radial flow. 
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FIG. 4: Transverse momentum spectra of dileptons for 
different invariant mass windows for minimum bias Au- 
Au collisions at RHIC energy. 
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FIG. 5: Variation thermal photon to dilepton ratio. Rem 
with Pt for different invariant mass windows at SPS en- 
ergy (see text). 



tion of the effective slope parameter extracted from 
the Mt spectra of the lepton pair with (M) eval- 
uated within the ambit of the present model [5^ 
reproduces the NA60 JET] results reasonably well. 
For Au-|-Au collisions at ^/snn=200 GeV, we have 
evaluated the dilepton spectra for different invariant 
mass bins with the initial condition (min bias) shown 
in Table I and lattice QCD equation of state. The 
results are displayed in Fig. S) The slopes of the ex- 
perimental data on px distribution of lepton pairs for 
different invariant mass windows measured by the 



B. The ratio, Ram 

As mentioned before some of the uncertainties pre- 
vailing in the individual spectra may be removed by 
taking the ratio. Rem of the px distribution of ther- 
mal photon to dileptons. In the absence of experi- 
mental data for both photon and dilepton from the 
same colliding system for SPS energies, we have cal- 
culated the ratio Rem for Pb-f Pb system, where the 
initial condition, the freeze-out condition and the 
EoS are constrained by the measured WA98 photon 
spectra. The results are displayed in Fig. [H 
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FIG. 6: Variation thermal photon to dilepton ratio, R^m 
with Pt for different invariant mass windows at RHIC 
energy (see text). 



FIG. 8: The variation of radial flow velocity with average 
temperature of the system for< M > =0.625 GeV and 
0.9 GeV at RHIC energy. 



Next we evaluate the ratio of the thermal photon 
to dilepton spectra constrained by the experimental 
data from Au+Au collisions measured by PHENIX 
collaboration. The results for the thermal ratio, Rem 
displayed in Fig. [6] is constrained by the experimen- 
tal data. The behaviour of i?om with pT for different 
invariant mass windows which is extracted from the 
available data is similar to the theoretical results ob- 
tained in Ref. The ratio, Rem for different M 
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FIG. 7: The variation of radial flow velocity with average 
temperature of the system for< M >= 0.75 GeV and 1.2 
GeV at SPS energy. 



windows (Figs. [5] and [6]) can be parametrized as 



follows: 



R, 



^ (MtY (Mt 

em = A exp -— 



PT 
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with different values of Teff for different invariant 
mass windows. The argument of the exponential in 
Eq. Hcan be written as [ll|; 



Mt 



Mt 
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Teff2 Teffl 
Mt 



Pt 
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1+tV 



(5) 



is the blue shifted effec- 



where, T^ffi = Ta 

five temperature for massless photons and Teff2 = 
Tav + Mv"^, is the effective temperature for massive 
dileptons. Tav is the average temperature and Vj. is 
the average radial flow of the system. For a given 
Pt and M Eq.[S]can be written as Vr — f{Tav)- The 
Tf,ff obtained from the parametrization of ratio at 
SPS energy are 263 MeV and 243 MeV for M=0.75 
and 1.2 GeV respectively. The average flow velocity 
Vr versus Tav have been displayed for M=0.75 GeV 
and 1.2 GeV in Fig.[ll The hadronic matter (QGP) 
dominates the M ~ 0.75(1.2) GeV region. There- 
fore, these two mass windows are selected to extract 
the flow parameters for the respective phases. The 
Vr increases with decreasing Tai, (increase in time) 
and reaches its maximum when the temperature of 
the system is minimum, i.e., when the system at- 
tains T/, the freeze-out temperature. Therefore, the 
variation of Vr with may be treated as to show 
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how the flow develops in the system. The Vr is larger 
in the hadronic phase because the velocity of sound 
in this phase is smaller, which makes the expansion 
slower as a consequence system lives longer - allow- 
ing the flow to fully develop. On the other hand, Vr 
is smaller in the QGP phase because it has smaller 
life time where the flow is only partially developed. 
In Fig. [5] the variation of average transverse veloc- 
ity with average temperature for RHIC initial condi- 
tions is depicted. The magnitude of the flow is larger 
in case of RHIC than SPS because of the higher ini- 
tial pressure. Because of the larger initial pressure 
and QGP life time the radial velocity for QGP at 
RHIC is larger compared to SPS. 

Obtaining Tc//! and Tcff2 from the individual 
spectra and eliminating Tat, one gets the variation 
of Vr with M. Fig. IH] (left panel) shows the varia- 
tion of Vr with M for SPS conditions. The radial 
flow velocity increases with invariant mass M up to 
M = Mp then drops. How can we understand this 
behaviour? From the invariant mass spectra it is 
well known that the low M (below p mass) and high 
M (above (j) peak) pairs originate from a partonic 
source [13] ■ The collectivity (or flow) does not de- 
velop fully in the QGP because of the small life time 
of this phase. Which means that the radial velocity 
in QGP will be smaller for both low and high M . 
Whereas the lepton pairs with mass around p-peak 
mainly originate from a hadronic source (at a late 
stage of the evolution of system) are largely affected 
by the flow resulting in higher values of flow velocity. 
In summary, the value of Vr for M below and above 
the p-peak is small but around the p peak is large 
- with the resulting behaviour displayed in Fig. [31 
Similar non-monotonic behaviour is observed in case 
of elliptic flow of photon as a function of pr (131 • The 
variation of Vr with M in RHIC (Fig. |9] right panel) 
is similar to SPS, though the values of Vr at RHIC 
is larger than that of SPS as expected due to higher 
initial pressure. 



laborations have been studied. The initial condi- 
tions have been constrained to reproduce the mea- 
sured multiplicity in these collisions. The EoS, the 
other crucial input to the calculations has been taken 
from lattice QCD calculations. The deviation of the 
hadronic phase from chemical equilibrium is taken 
in to account by introducing non-zero chemical po- 
tential for each hadronic species. It is shown that 
simultaneous measurements of photon and dilepton 
spectra in heavy ion collisions will enable us to quan- 
tify the evolution of the average radial flow velocity 
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FIG. 9: The variation of radial flow with invariant mass 
pairs for SPS (left) and RHIC (right) energies. 



for the system and the nature of the variation of ra- 
dial flow with invariant mass indicate the formation 
of partonic phase at SPS and RHIC energy. 



V. SUMMARY AND CONCLUSIONS 

The photon and dilepton spectra measured at SPS 
and RHIC energies by different experimental col- 
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